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^ . Recent neutrino data have been favourable to a nearly bimaximal mixing, 



which suggests a simple form of the neutrino mass matrix. Stimulated by this 

matrix form, a possibility that all the mass matrices of quarks and leptons 

have the same form as in the neutrinos is investigated. The mass matrix 

form is constrained by a discrete symmetry Z3 and a permutation symmetry 

S2. The model, of course, leads to a nearly bimaximal mixing for the lepton 

sectors, while, for the quark sectors, it can lead to reasonable values of the 

CKM mixing matrix and masses. 
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I. INTRODUCTION 



Recent neutrino oscillation experiments [1] have highly suggested a nearly bimaximal 
mixing (sin^ 2^i2 ~ 1, sin^2^23 — 1) together with a small ratio R = Amf2/Am|3 ~ 10~^. 
This can be explained by assuming a neutrino mass matrix form [2]- [7] with a permutation 
symmetry between second and third generations. We think that quarks and leptons should 
be unified. It is therefore interesting to investigate a possibihty that all the mass matrices 
of the quarks and leptons have the same matrix form, which leads to a nearly bimaximal 
mixing and C/13 = in the neutrino sector, against the conventional picture that the mass 
matrix forms in the quark sectors will take somewhat different structures from those in the 
lepton sectors. In the present paper, we will assume that the mass matrix form is invariant 
under a discrete symmetry Z3 and a permutation symmetry S2. 

Phenomenologically, our mass matrices M„, M^, My and Mg (mass matrices of up quarks 
(ti, c, t), down quarks {d,s,b), neutrinos {i'c^hj^t) and charged leptons {e,iJ,,T), respec- 
tively) are given as follows: 

Mf = PlfMfPnf, (1.1) 

with 



Mf 



^ Af Af^ 

Bf Cf 
\ Cf Bf J 



{f^u,d,u,e), (1.2) 



where Pl/ and Pr/ are the diagonal phase matrices and Af, Bf, and Cf are real parameters. 
Namely the components are different in Mf, but their mutual relations are the same. This 
structure of mass matrix was previously suggested and used for the neutrino mass matrix in 
Refs [2]- [7], using the basis where the charged-lepton mass matrix is diagonal, motivated 
by the experimental finding of maximal i^^^-i^r mixing [1]. In this paper, we consider that 
this structure is fundamental for both quarks and leptons, although it was speculated from 
the neutrino sector. Therefore, we assume that all the mass matrices have this structure. 



Let us look at the universal characters of the model. Hereafter, for brevity, we will omit 
the flavour index. The eigen-masses mj of Eq. (1.2) are given by 



-mi = ^ + C - ^J^¥V{BTcf^ , 
m2 = ^ + C + ^J^¥V{B^^^ , 



The texture's components of M are expressed in terms of eigen-masses rrii as 



A 
B 



1 ^2 - mi 
-ma 1 + 

2 \ ms 



^ 1 / m2 - mi 

C= --ma 1 

2 V ma 



That is, M is diagonalized by an orthogonal matrix O as 

^-mi 0^ 
0^M0= m2 

^0 1713 J 



with 



O = 



c s 

s_ _c_ 1_ 

s_ _c_ 1 

\ V2 V2 V2 / 



Here c and s are defined by 



m2 



mi 



(1.3) 
(1.4) 
(1.5) 



(1.6) 



[1.7) 



(1.8) 



(1.9) 



y m2 + mi y m2 + mi 

It should be noted that the elements of O are independent of ms because of the above 
structure of M. 

The zeros in this mass matrix are constrained by the discrete symmetry that is discussed 
in the next section, defined at a unification scale (the scale does not always mean "grand 



unification scale"). This discrete symmetry is broken below fx — Mr, at which the right- 
handed neutrinos acquire heavy Majorana masses, as we discuss in Sec. IV. Therefore, the 
matrix form (1.1) will, in general, be changed by renormalization group equation (RGE) 
effects. Nevertheless, we would like to emphasize that we can use the expression (1.1) with 
(1.2) for the predictions of the physical quantities in the low-energy region. This will be 
discussed in the appendix. 

This article is organized as follows. In Sec. II we discuss the symmetry property of our 
model. Our model is realized when we consider two Higgs doublets in each up-type and 
down-type quark (lepton) mass matrices. The quark mixing matrix in the present model is 
argued in Sec. III. In Sec. IV, the lepton mixing matrix is analyzed. Sec. V is devoted to a 
summary. 



II. Zs SYMMETRY AND MASS MATRIX FORM 



We assume a permutation symmetry between second and third generations, except for the 
phase factors. However, the condition (M/)ii = cannot be derived from such a symmetry. 
Therefore, in addition to the 2 3 symmetry, we assume a discrete symmetry Z3, under 
which symmetry the quark and lepton fields ipL, which belong to IOl, 5l and 1l of SU(5) 
{1l — J^r), are transformed as 

where cu^ — +1. (Although we use a terminology of SU(5), at present, we do not consider 
the SU(5) grand unification.) Then, the bilinear terms qj^^URj, q^^dRj, iu^^Rj, iu^Rj and 



^Ri^Rj [^R = {^rY = Cur^ and u'^r = (i^r)] are transformed as follows: 

1 UJ'^ OJ^ 



LO LO LO 



LO LO LO 



(2.2) 



where 



Ql = 



(2.3) 



Therefore, if we assume two SU(2) doublet Higgs scalars Hi and H21 which are transformed 
as 



Hi — > ujHi^ H2 — > Lo'^H2, 
the Yukawa interactions are given as follows 

-f^int = {^{A)ij^LiHAURj + Y^^j^^-j-q^iHAdRj) 

A=l,2 

+ X! (y(A)i/LiHA^Rj + Y^^yjIuHAeRj^ 
A=l,2 



(2.4) 



(2.5) 



where 



Ha 





, Ha = 


( If \ 

^A 








[-Ha) 



(2.6) 



so that 



^0 0^ 



V" V" — 
^(2), ^(1), ^(2), -^(2) - 



* * 
* * 



u * * 



VU -i^e \^R 

^(2), ^{ly ^(2), ^(1), ^(1) 



* 

* 



(2.7) 



In (2.7), the symbol * denotes non-zero quantities. Here, in order to give heavy Majorana 
masses of the right-handed neutrinos vr, we have assumed an SU(2) singlet Higgs scalar 
which is transformed as Hi. 

In the present model, the phase difference arg(Y"^{^ -|- Y^^-^)2i — a'rg(Y^{^ -|- y^2))^'^ plays 
an essential role. Therefore, for the permutation symmetry S2, we put the following as- 
sumption: the permutation symmetry can be applied to only the special basis that the all 
Yukawa coupling constants are real. (Of course, for the Z3 symmetry, such an assumption 



is not required.) We consider that the phase factors are caused by an additional mechanism 
after the requirement of the permutation symmetry S2 (after the manifestation of the hnear 
combination + Y(2)). In the present paper, we consider that although the Z3 symmetry is 
rigorously defined for the fields by (2.1), the permutation symmetry S2 is rather phenomeno- 
logical one (i.e. Ansatz) for the mass matrix shape. Then, under such the S2 symmetry, the 
general forms of Yf = Y,L + Y,L are given by 



(1) -r - (2) 

Qf>—i{S{^2~^Rl) l)(^~^i^L2~^R2) Ce~*('^i2~'^i{3) 



y ae-'^^is-^Ri^ ce-*('^M-42) fte^^^'^M-'^fls) j 



(2.8) 



We have already assumed that ipL — {j^,eL,d%;uL,dL,u'^ji,e'^ji; J^r) have the same transfor- 
mation (2.1) under the discrete symmetry Z3, so that ipLii'lj are transformed as (2.2). Prom 
this analogy, we assume that the phase matrices P^f and Pjif come from the replacement 
tpL Pfi^L, i-e. 

^^Yf^L ^ ^LPjYfPjrL ■ (2.9) 

However, differently from the transformation (2.1), we do not assume in (2.9) that all the 
phase matrices Pf are identical, but we assume that they are flavour dependent. This 
explains the assumption 

SL = -4 = 5/ , (2.10) 

in the expression (2.8). (However, this assumption (2.10) is not essential for the numerical 
predictions in the present paper, because the predictions of the physical quantities depend 
on only the phases (5{j.) 

Since the present model has two Higgs doublets horizontally, in general, flavour- changing 
neutral currents (FCNCs) are caused by the exchange of Higgs scalars. However, this FCNC 
problem is a common subject to be overcome not only in the present model but also in most 
models with two Higgs doublets. The conventional mass matrix models based on a GUT 
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scenario cannot give realistic mass matrices without assuming more than two Higgs scalars 
[8]. Besides, if we admit that two such scalars remain until the low energy scale, the well- 
known beautiful coincidence of the gauge coupling constants at ~ 10^^ GeV will be spoiled. 
Although the present model is not based on a GUT scenario, as are the conventional mass 
matrix models, for the FCNC problem, we optimistically consider that only one component of 
the linear combinations among those Higgs scalars survives at the low energy scale = mz, 
while the other component is decoupled at < Mx [9] . The study of the RGE effects given 
in the appendix will be based on such an "effective" one-Higgs scalar scenario. 



III. QUARK MIXING MATRIX 



The quark mass matrices 



Mf = PjMfPj {f = u,d), 



(3.1) 



are diagonalized by the bi-unitary transformation 



Df = UlfMfURf , 



(3.2) 



where Uif = P}Of, Ur} = PfOf, and Od (0„) is given by Eq. (1.8). Then, the Cabibbo- 
Kobayashi-Maskawa (CKM) [10] quark mixing matrix V is given by 

V = UlULd = O^P^P^Od 



CuCd + pSuSd CuSd - pSuCd 

SuCd - pCuSd S^Sd + pCuCd 

-asd crcd 
where p and a are defined by 



ac, 



u 



V 



p=l(e'53 + e^5^) = cos^^ 



expz 



'^3 + 62^ 



(3.3) 



(3.4) 



2^ ^ 2 



expW^- + - 



(3.5) 



Here we have put P = -Pu-Pj = diag(e*''% e*'^^, e^^^), and we have taken 5i — without loss of 
generahty. 

Then, the exphcit magnitudes of the components of V are expressed as 



IK 



sm 



S3—S2 



cb 



a c,. 



IK 



y^l + m„/mc 
sin^^ 



ub\ 



\a\ a 



\Vts\ = |(7| Cd = 

\Vtd\ = |(j| = 



sm 



$3—52 



sin^^ 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



X 



1 o - (^2 ^3 + 
1 — Z COS COS 



mc + mu\l ms + rud 
3 



rricmd 



+ COS 



2 «3 



m„m. 



rricmd 



(3.10) 





Su Cd 






p 











rrid 



mc + muy ms + rrid 



X 



2 (^3 - (^2 „ 5^-82 5^ + §2 Frrrjus 
cos — 2 cos — - — cos — - — 4 / h 



(3.11) 



2 2 2 y rricmd \mcmd. 

It should be noted that the elements of V are independent of and m^. The independent 
parameters in the expression |Kj| arc 9u = tan~^(m„/mc), 9d = tan~^{md/ms), 63, and 62- 
Among them, the two parameters 9u and 9d are already fixed by the quark masses of the 
first and second generations. Therefore, the present model has two adjustable parameters 
^3 and 62 to reproduce the observed CKM matrix parameters [11]: 



IK 



us I cxp 



0.2196 ± 0.0026, iK^lexp = 0.0412 ± 0.0020, 



|K6|,.p = (3.6±0.7) X 10-^ 
It should be noted that the predictions 



(3.12) 



\Vcb\ 

\Vtd\ 
IK 



ts 



Sd 
Cd 



jmu 

mc 

ms 



/2.33 
677 

/4.69 
93A 



0.0586 ± 0.0064 , 



0.224 ±0.014 



(3.13) 
(3.14) 



are almost independent of the RGE effects, because tfiey do not contain the phase difference, 
(^3 — S2), which is highly dependent on the energy scale as we discuss in the appendix [see 
(A. 9)] and we know that the ratios m„/mc and md/mg are almost independent of the RGE 
effects. In the numerical results of (3.13) and (3.14), we have used the running quark mass 
at II — mz [12]: 



The predicted value (3.13) is somewhat small with respect to the present experimental value 
iK&l/lKbl = 0.08 ± 0.02, but it is within the error. 

The heavy-quark-mass-independent predictions (3.13) and (3.14) have first been derived 
from a special ansatz for quark mixings by Branco and Lavoura [13], and later, a similar 
formulation has also be given by Fritzsch and Xing [14]. For example, the CKM matrix V is 
given by the form V — Ri2{0u)R2z{&Q, (t>Q)F^2{^d) in the Fritzsch-Xing ansatz, and their ro- 
tation R2s{0q, 4>q) with a phase (j)Q corresponds to i?23(— 7'"/4)-Pu-Pj-R^3(— 7r/4) in the present 
model, because the present rotation given in (1.8) is expressed as 0/ = i?23(— 7r/4)i?i2(^/). 
However, we would like to emphasize that the 2 <-> 3 mixing in V comes from only the rela- 
tive phase difference {52 — 5^), and it is independent of the forms of the up- and down-mixing 
matrices (1.8). The present mass matrix texture is completely different from theirs. The 
rederivation of (3.13) and (3.14) in the present model will illuminate the farsighted instates 
by Branco and Lavoura. 

Next let us fix the parameters ^3 and 52- When we use the expressions (3.6)-(3.11) at 
IX — mzi the parameters 82 and ^3 do not mean the phases that are evolved from those 
Ski II — Mx- Hereafter, we use the parameters 82 and ^3 as phenomenological parameters 
that approximately satisfy the relations (3.6)-(3.11) at = mz- In order to fix the value of 
5^ — 82, we use the relation (3.6), which leads to 



mu{mz) = 2.33lHiMeV, mc{mz) = 677l|?MeV, 
md{mz) = 4.69lg|^MeV, m,(mz) = 93.4ll^:« MeV. 



(3.15) 




(3.16) 



4.59° ± 0.21°. 



(3.17) 
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Then, we obtain 



iKbl 



im. 



rrir 



cb\ 



1 + 



\1 + 



— \Vcb\ 
md I ™ I exp 



0.00234 ± 0.00028, 
0.0391 ± 0.0018, 



w. 



td\ 



1 _i_ m .. 



-<—Kb\ = 0.00880 ± 0.00094, 



rus' ^""""^ 



(3.18) 
(3.19) 
(3.20) 



which are consistent with the present experimental data. Therefore, the value (3.17) is 
acceptable as reasonable. Then, by using the value (3.17) and the expression (3.10), we can 
obtain the remaining parameter (^3 + ^2) : 



53 + 82^ 93° ± 22° or - 80° ± 22° . 



(3.21) 



Since sin(53 — 62) /2 ~ 0.04 and cos((53 + S2)/2 ~ 0.2, the present model also predicts the 
following approximated relations 



IK 

\Vcd\ = CuSd 



l-p 



Irud 



Cd 

p 

Cti Sd 



IrUd 



\Vtd\ = \cr\ Sd = \/\Vcbf + \ Vub\ 



rrid 



Wcbl ■ \Vus\ 



(3.22) 
(3.23) 
(3.24) 



nis + md 

Using the rephasing of the up-type and down-type quarks, Eq. (3.3) is changed to the 
standard representation of the CKM quark mixing matrix 



^4 . 



Ktd = diag(e"i , ,e^^)V diag(e"i , e"^ , e"^ ) 



C13C12 



C13S12 



si3e 



V 



-C23S12 - S23Cl2Sl3e C23C12 - S23Sl2Sl3e S23C13 



(3.25) 



S23S12 - C23Ci2Si3e -S23C12 - C23Sl2Si3e C23C13 

Here, af comes from the rephasing in the quark fields to make the choice of phase convention. 
The CP-violating phase 5 in the representation (3.25) is expressed with the expression V in 
Eq. (3.3) by 
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S — arg 

so that we obtain 





Vl2 


2 ' 


1 - 






2 



(3.26) 



5 = ±(80° ±22°). (3.27) 

It is interesting that nearly maximal | sin 5\ is realized in the present model. 

The rephasing invariant Jarlskog parameter J [15] is defined by J = Im{VusV*gV*jjVcb) ■ 
In the present model with (3.6)-(3.11), the parameter J is given by 

T I |2| I . ^3+^2 

J=|(7| \p\CuSuCdSdSm — - — 

sm — - — . (3.28) 

Using the relation \ytd\ — iKbllKsl in (3.24), and the experimental findings |T4sP ^ iKbP ^ 
iKfcP, \Vts\ ^ \Vcb\, and \Vtb\ ~ 1, we obtain 

J~ iKfellV^fellK^lsin^^^. (3.29) 

On the other hand, in the standard expression of V, (3.25), J is given by 

J = C^3Si3Ci2Si2C23S23 siu S 



\Vub 




Vtd 


\Vts\ 


\Vtb\ 


\Vcb 


l+\Vub/Vcb\^ 



-13" 

\Vud\ IKjsI iKjbl iKbl l^tbl 
12 



sin 5 ~ I Vus 1 1 Kb 1 1 K& I sin 5. (3. 30) 



l-\Vub\ 

Comparing Eq. (3.29) with Eq. (3.30), we obtain 

. r . S3 + S2 , . 

sm()~sm — - — . (3.31) 

By using the numerical results (3.17)-(3.21), we obtain 

|J| = (1.91±0.38) X 10-^ (3.32) 
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IV. LEPTON MIXING MATRIX 



Let us discuss the lepton sectors. We assume that the neutrino masses are generated via 
the seesaw mechanism [16]: 



(4.1) 



Here Md and Mr are the Dirac neutrino and the right-handed Majorana neutrino mass 
matrices, which are defined by VlM^vr and V^Mri/r, respectively. Since Mo — P^MoP^ 
and Mr — P^MrP}, according to the assumption (2.9), we obtain 

= -pImdM^^mIpI 

I 



Pi 







m2m 



I m2mi [ 

V 2 V 

>3(l + «) -\m,{l 

^->3(l-^) >3(l- 



m2mi 
2 



TO2mi 
2 2 



m2— wii 

m2— mi 
ma 



Pi . 



(4.2) 



2 2"'''^V rns J 2' 

Here and hereafter, mi, m2 and 777,3 denote neutrino masses unless they are specifically 
mentioned. In the last expression, we have used the fact^ that the product of AB~^A of the 
matrices A and B with the texture (1.1) with (1.2) again becomes a matrix with the texture 
(1.1) with (1.2). 

On the other hand, the charged lepton mass matrix Me is given by 



Me = Pi 







m^me 1 



mr{l + ^^^^ 



-2^r [I - 







e ' 



(4.3) 



where me, m^ and m,- are charged lepton masses. 

Those mass matrices M^ and My are diagonalized as {PlOe^ M^^PeOe 
(Pj;0,,)^M^(Pj,Oj,) = i:)^, respectively, where 



Df, and 



The seesaw invariant texture form was discussed systematically in [17]. 
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Ce 


Se 





Oe = 




Ce 


1 


V2 


V2 


V2 




Se 

\ V2 


Ce 


1 




V2 


%/2 














1 


V2 




V2 


s„ 




1 


V2 


V2 


a/2 



(4.4) 



Here Cg and Sg are obtained from Eq. (1.9) by replacing mi and m2 in it by rUe and m^; c,^ 
and Sn are also obtained by taking the neutrino masses mj. Therefore, the Maki-Nakagawa- 
Sakata-Pontecorv (MNSP) lepton mixing matrix [18] U can be written as 



U = Oi PO^ 
I 



SeCu — pyCeSu SeS^ + p 
-OuS 



OvCe 
Pv 



(4.5) 



where P = PePl = diag(e*''''i, e^'^''^, e*''''^). Hereafter we take 5i,i — without loss of general- 
ity. 

The explicit forms of absolute magnitudes of the components of U are given by expres- 
sions similar to (3.4)-(3.12), where \Vij\, {mu, rric, rrit), and (m^, m^, nib) are replaced by \Uij\, 
{mi, 1712,1713), and {me,mn,mr), respectively. It should again be noted that the elements 
of U are independent of rrir and 7713. The independent parameters of the unitary matrix 
U are 6^ = tan~^ (mg/m^), 61, — tan~^(mi/m2), 8^3, and 5^2- Among them, 6^ is given by 
charged-lepton masses of the first and second generations. Therefore, the model has the 
three adjustable parameters 5i,z, 5^2, and mi/ 1712 to reproduce the experimental values [11]. 

Let us estimate the values 9i,, 5i,3 and 5,^2 by fitting the experimental data. In the 
following discussions we consider the normal mass hierarchy Am^a = m| — ml > for the 
neutrino mass. The case of the inverse mass hierarchy Am^g < is quite similar to it. It 
follows from the CHOOZ [19], solar [20], and atmospheric neutrino experiments [1] that 



Prom the global analysis of the SNO solar neutrino experiment [20] , 



(4.6) 
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Am?2 ^ml-ml^ Aml^ = 5.0 x 10~^ eV^ (4.7) 

tan^ = tan^ 6,,^ = 0.34, (4.8) 

with x^j^/dof = 57.0/72, for the large mixing angle (LMA) MSW solution. From the 
atmospheric neutrino experiment [1], we also have 

Am^3 = - ~ Am^^^ = 2.5 x 10"^ eV^ (4.9) 

sin" 2^23 ^sin^ 2^,^^ = 1.0, (4.10) 

with xLn/dof = 163.2/170. 

Independently of the parameters 6^3 and 61,2, the model predicts the following two ratios: 



If^isl 


Se 1 


\U23\ 




\U3l\ 


Su / 


\U32\ 





0.0688, (4.11) 



V 103 



^ . (4.12) 

m2 

Here we have used the running charged-lepton mass at /i = mz [12]: meijnz) = 0.48684727± 
0.00000014 MeV, and m^{mz) = 102.75138 ± 0.00033 MeV. The neutrino mixing angle 6,,^ 
under the constraint [Amlgl ^ |Amf2| is given by 

sin2 2^,,^ = 4|C/23|'|C/33|' 



= 4|p.|2 \a,fcl = sm\6,s-6,2)J ■ (4.13) 

^ m^ + me 

The observed fact sin^2^atm — 1.0 highly suggests 6^3 — — 7r/2. Hereafter, for simphcity, 
we take 

S.3 -^u2 = \. (4.14) 
Under the constraint (4.13), the model predicts 

1 Tfl 

l^/ial^ = — = 0.00236 , or sin22^i3 = 0.00942 . (4.15) 

2 + me 

This value is consistent with the present experimental constraints (4.6) and can be checked 
in neutrino factories [21], which have sensitivity to sin^2^i3 for 
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8111^2^13 > 10-^ 



(4.16) 



The mixing angle 9^„i in the present model is given by 



sin2 2^.„, = 4|C/nnC/i2| 
4m2mi 



which leads to 



(m2 + mi)^ 
4mi/m2 
(1 + mi/m2Y 



\-\f2 



cos ■ 



+ 



1 f mf.m2\ 



— ~ tan^ ^3oi = 0.34 , 

777.2 



(4.17) 



(4.18) 



where we have used the best fit value (4.8). This value (4.18) guarantees the vahdity of 



the approximation (4.17), because of \J {rrie/ m^) / {mi/ — 0.12. Then, we can obtain the 
neutrino masses 



mi = 0.0026 eV , 
m2 = 0.0075 eV , 
ma = 0.050 eV , 



(4.19) 



where we have used the observed best fit values of Am^^j and Am^^.^, (4.7) and (4.9), respec- 
tively. 

Next let us discuss the CP violation phases in the lepton mixing matrix. The Majorana 
neutrino fields do not have the freedom of rephasing invariance, so that we can use only the 
rephasing freedom of Mg to transform Eq. (4.5) to the standard form 



Ustd = 



JI3 



Cu23Sul2 - S:.23C:.23Si/13e"^")e~'^ C^23C^12 " S,.23S,.12S,.13e"^" 8^230^13^'^'^''^^ 

-Su23Cvl2 - Cu23Si,12Sul3e'^'')C^^~^^ Cu23Cul3 



iSu23Sul2 - Cu23Cul2Sul3e'^")e ^'^ 



(4.20) 



as 
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Ustd = diag(e^"i, e^"2, e'"^) U diag(e±^^/^ 1, 1) . (4.21) 

Here, a? comes from the rephasing in the charged lepton fields to make the choice of phase 
convention, and the specific phase ±7r/2 is added on the right-hand side of U in order to 
change the neutrino eigen-mass mi to a positive quantity. Similarly to the quark sector, the 
CP-violating phase in the representation (4.20) is expressed as 



5u = arg 



U12U22 _|_ 





u 


12 


2 


1 - 




[/l3 


2 



-argi ^y-^ 1 ~argp, + 7r = ^ (4-22) 



.[/l3[/|3 

Though the lepton mixing matrix includes the additional Major ana phase factors j3 and 
7 [22,23], the number of parameters which will become experimentally available in the 
near future is practically four, as in the Dirac case. The additional phase parameters are 
determined as 

and 

by using the relations me <S m^ and {5,^3 — Si,2)/2 ~ 7r/4. Hence, we can also predict the 
averaged neutrino mass {m,y) [23], which appears in the neutrinoless double beta decay, as 
follows: 

(m^) = -miUli + m2t/i2 + "^3t^i3 

= -2p^CeS^y/rnrm^ + plsl{m2 - rui) +mssl . (4.25) 

The value of (4.25) is highly sensitive to the value of (8^3 + 61,2) f^, which is unknown at 
present, because the values Se/ce = \Jrne/m^ ~ 0.070 and y^mTm^/ms ^ 0.088 are in the 
same order. For {5„3 -\- 6„2)/'2 — 0, 7r/2 and tt, we obtain the numerical results {rrii) — 
0.00018 eV, 0.00049 eV and 0.00069 eV, respectively. However, these values should not be 
taken strictly because the value mi/m2 is also sensitive to the observed value of tan^^so^. 
In any cases, the predicted value of (m,^) will be less than the order of 10~^ eV. 
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The rephasing-invariant parameter J in the lepton sector is defined hy J — 
Im{Ui2U22U^3U23) , which is exphcitly given by 

J = \<^u\ \Pu\CuS,yCeSeSm ■ 



\Ul3\ 


\U3l\ 


|f/32| 




\U23\ 


l + lf/isMsP 



^1/3 ~l~ ^1^2 , I ^ ici I I ^ I I ^ ,5Z I I ^ I /I r>/-\ 

sm < -J- — ' ' . (4-26) 



\Ul3\ 




\U3l\ 


\U32\ 


\U33\ 


\U23\ 


1 + 1 


U13/U23V 



2 

3_ 

~2 

The upper bound is described in terms of the ratio mi/m2, so that we obtain 

J < 0.019. (4.27) 

It should be noted that if we again assume the maximal CP violation in the lepton sector, 
the magnitude of the rephasing invariant \J\ can be considerably larger than in the quark 

sector, I JquarkI — 2 X 10~^. 



V. CONCLUSION 

In conclusion, stimulated by recent neutrino data, which suggest a nearly bimaximal 
mixing, we have investigated a possibility that all the mass matrices of quarks and leptons 
have the same texture as the neutrino mass matrix. We have assumed that the mass matrix 
form is constrained by a discrete symmetry Z3 and a permutation symmetry S2, i.e. that 
the texture is given by the form (1.1) with (1.2). The most important feature of the present 
model is that the textures (1.1)-(1.2) are practically applicable to the predictions at the low 
energy scale (the electroweak scale) , although we assume that the textures are exactly given 
at a unification scale. 

It is well known that the matrix form (1.1) leads to a bimaximal mixing in the neutrino 
sector. In the present model, the mixing angle ^{2 between the first and second generations 
is given by 

tan^{2 = ^(A4, (5.1) 

where ra{ and ra^ are the first and second generation fermion masses. This leads to a 
large mixing in the lepton mixing matrix (MNSP matrix) U with mi ~ m2 (neglecting 
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tan^f2 — Y^?7T.e/m^ in the charge lepton sector), and it also leads to the famous formula 



[24] \Vus\ — \Jfnci/ms in the quark mixing matrix (CKM matrix) V (neglecting tan^^"2 = 
niu/ rric in the up-quark sector). In the present model the mixing angle ^4 between the 
second and third generation is fixed as ^4 — However, the (2,3) component of the 

quark mixing matrix V (and also the lepton mixing matrix U) is highly dependent on the 
phase difference S3 — 62, as follows 

V23 = y sm -—— , (5.2) 

Jl + mi/m2 ^ 

where 5j = 5" — 5f . Replacing the arguments by their leptonic counterparts, we have the 
same form for C/23- We have understood the observed values V23 and [/23 by taking ib^ — b-i) /2 
as a small value for the quark sectors and as 7r/2 for the lepton sectors, respectively. As 
predictions, which are independent of such phase parameters, there are two relations 



IKbl \yu\ 




(5.3) 



(and the similar relations for [/) . The relations (5.3) are in good agreement with experiments. 



The relation IC/13/C/23I = yme/m^ in the lepton sectors leads to ^ me/2m^ — 0.0024 

if we accept sin^26'atm = 1-0. This value will be testable in the near future. 

Since, in the present model, each mass matrix Mf (i.e. the Yukawa coupling Yf) takes 
different values of Af, Bf, and so on, the present model cannot be embedded into a GUT 
scenario. In spite of such a demerit, however, it is worth while noting that it can give a 
unified description of quark and lepton mass matrices with the same texture. 
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The mass matrix texture (1.1) with (1.2), which is defined at the unification energy scale 
IJ, = Mx, is applicable to the phenomenology at the electroweak scale = mz- In the 
present appendix, we demonstrate this for the quark mass matrices and M^. 

It is well known [25] that the energy scale dependences R{A) — A{ijl)/A{Mx) for ob- 
servable quantities A approximately satisfy the relations R{\Vub\) — R{\Vcb\) — R{\Vtd\) — 
R{\Vts\) — R{md/mb) — R{ms/mb), and that the ratios i?(|V^s|), i?(|V^d|), R{jnd/ms) and 
R{mu/m^ are approximately constant. This is caused by the fact that the Yukawa coupling 
constant yt of the top quark is extremely large with respect to other coupling constants. The 
above relations on R are well explained by the approximation y| ^ yl-iV^i' ' Therefore, 
we will also use approximation below. 

The one-loop RGE for the Yukawa couphng constants 1/ (/ = u, d) has the form 

dYf _ 



{Cfl + C,,Y,Y} + C,,Y,Y}) Y, , 



(Al) 



dt 167r2 

where f' — d (/' — u) ior f — u {f — d), and the coefficients C/, C// and C/// are energy 
scale dependent factors which are calculated from the one-loop Feynman diagrams. Wc start 
from the Yukawa coupling constants l/(Mx), corresponding to the mass matrix form (1.1) 
[with (1.2)]. 

Since the matrix YyY^ is approximately given by 

/ 



Y^{Mx)Yl{Mx) 




1 
-e-' 



— e 







\ 



{A.2) 



I 



where 5„ = ^3 — ^2, Vuj = {H^), and we have used the relations (1.6) and the approxima- 
tion Ut ^ ybTVcT ' ' '■> up-quark Yukawa coupling constant l^(/x) in the neighbourhood of 
II — Mx is given by the form 



y„(//) ~ r„(//) 





1 -e+'^^ 
-e-'^^ 1 



Yu{Mx) 
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A^e-'^'^ A^e-^^l 



(A3) 

Although this form is one in ~ Mx, but, since the texture keeps the same form under 
the small change of energy scale, as a result, the texture of given by (A. 3) holds at 

any energy scale /x. Therefore, we can obtain the expression (1.1) at an arbitrary energy 
scale \x. (The demonstration (A. 3) has been done for the case Pr = mentioned in (2.10). 
However, the conclusion does not depend on this choice.) 

On the other hand, the evolution of the down-quark Yukawa coupling constant Yd{,y^ is 
somewhat complicated. By a way similar to (A. 3), we obtain 











l+£d 



^ -e^e-'^- 1 + £d y 



( 



Aa 
Aa{l + £d - £,e+^('^-*'')) Bail + £d - Sae+'^'^-^'^^Ca/ B^) C,(l + - £<ie+^('^— ^<^)S,/Crf) 



^ Ad{l + ed- ede-"^^'--'-y) Crf(l + - sae-'^'--''^^ Ba/Ca) 5^(1 + 6^- ede~'^'--'^)Cd/ B^) 

[AA) 

where 5u — 5d = {5^ — 5f) — (5| — Sf) = 5^ — 52- Note that the part that is sandwiched 
between P\ and P\ includes imaginary parts and those phase factors cannot be removed by 
an additional phase matrix Pd(/x) into the form Pj(/x)yrf(/i)Pj(yu). However, the quantity 
that has the physical meaning is Ydf^. When we define 



-i{Su-Sd) 



) , Ve 



1 + eJl + e 



+i{S„-Sd) 



{A.5) 



20 



we obtain 



42 



Ad{B, + C,)r) Aie + {Bj + Cl)rf A\e^^^-f^) + 2^,^,77' 



(A.6) 



where 



(A.7) 



so that we can obtain a real matrix for the part which is sandwiched by the phase matrix 
under the approximation Ai^J\BdC^ ~ 0. This means that we can practically write 



yM - 



Vd/V2 



Ad Ad 
AdC BdT] CdV 
AdC CdV BdV 



\ 



{A.8) 



with 



Pj(//) = diag(l,e-H^^"-^\e 



{A.9) 



at an arbitrary energy scale /i. It should be noted that the changes of the phases 62^62—^ 
and (5f — > 5f + /3 do not come from the evolution of the phases (^2 (a*) ^iip') 1 but they 
are brought effectively by absorbing the unfactorizable phase parts in Yrf(/x). Thus, we can 
again use the texture (1.1) at an arbitrary energy scale ^ from a practical point of view. 

In the Yukawa coupling constants Yg and Yj^ of the leptons, the RGE effects are not so 
large as in the quark sectors. In the charged lepton sector, since » ^ ml, we can 
again demonstrate that the expression (1.1) is applicable at an arbitrary energy scale in a 
way similar to the quark sectors. For the neutrino Yukawa coupling constant 1^(//), the 
evolution equation is different from (A.l). We must use the RGE for the seesaw operator 
[26]. However, the calculation and result are essentially the same as those in F„(/x), Yd(/x) 
and Ye{^), because m\ ^ m| > m{ in the present model. 

Finally, we would hke to add that these conclusions on the evolution of the mass matrices 
M/ (/ = M, (i, e, p) are exactly confirmed by numerical study, without approximation. 
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